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We present a fully nonlinear theory for dust acoustic (DA) shocks and DA solitary pulses in a 
strongly coupled dusty plasma, which have been recently observed experimentally by Heinrich et 
al. [Phys. Rev. Lett. 103, 115002 (2009)], Teng et al. [Phys. Rev. Lett. 103, 245005 (2009)], 
and Bandyopadhyay et al. [Phys. Rev. Lett. 101, 065006 (2008)]. For this purpose, we use 
a generalized hydrodynamic model for the strongly coupled dust grains, accounting for arbitrary 
large amplitude dust number density compressions and potential distributions associated with fully 
nonlinear nonstationary DA waves. Time-dependent numerical solutions of our nonlinear model 
compare favorably well with the recent experimental works (mentioned above) that have reported the 
formation of large amplitude non-stationary DA shocks and DA solitary pulses in low-temperature 
dusty plasma discharges. 

PACS numbers: 52.27.Lw,52.35.Tc,52.35.Fp 



I. INTRODUCTION 

Charged dust grains and dusty plasmas [IHZ] are ubiquitous in astrophysical environments (e.g. interstellar media, 
molecular dusty clouds, star forming clouds, supernovae such as the Eagle Nebula, etc.), in planetary ring systems [TJIH] 
(e.g. the spokes in Saturn's rings recorded by the Voyager spacecraft cameras), in our solar system (e.g. interplanetary 
dust particles produced by comets), as well as near the Sun's and Earth's atmospheres (e.g. the mesospheric and 
ionospheric regions). Charged dust particles are naturally formed in industrial processing of nanotechnology and in 
magnetic fusion reactors. 

It is well-known that charging of a neutral dust particle occurs due to a variety of physical processes [51 UOj , including 
the collection of electrons from the background plasma, photo emissions, tribo-electric effects, etc. In the remote past, 
it was shown by Wuerker et al. |11) that an ensemble of electrically charged iron and aluminum particles having 
diameters of a few microns can be confined by three-dimensional focusing forces of alternating and static electric 
fields and the Coulomb repulsion, leading eventually to the formation of crystallized arrays of ions and aluminum dust 
particles, which can be melted and reformed. However, a dusty plasma is usually composed of electrons, positive ions, 
negative or positive dust grains, and neutral atoms. When the interaction potential energy (— Z 2 e 2 /d, where Zd is 
the dust charge state, e the magnitude of the electron charge, and d the inter-dust grain distance or the Wigner-Seitz 
radius) between two neighboring dust grains is much larger (smaller) than the dust kinetic energy ksTd, where ks 
is the Boltzmann constant and Td the dust temperature, the dusty plasma is in a strongly (weakly) coupled state. 
Following the charged particles condensation idea (T2] of one component strongly correlated electron system, Ikezi [T3] 
postulated the solidification of charged dust particles when the dusty plasma T = Z 2 d e 2 exp(— <i/ 'Ad) / 'dksTd exceeds 
172, taking into account the plasma screening effect, where Ad is the plasma Debye radius [2]. Such values of T 
can be achieved in low-temperature laboratory discharges at room temperatures owing to the large Zd acquired by a 
micron-size dust grain by absorbing electrons from the background plasma. There are also Monte-Carlo and Molecular 
Dynamics simulations that accurately depict different states of ordered dust structures [TlrfTrj] when dust grains are 
repelling each other according to the Yukawa or Debye-Huckel force. The phase diagram for T s against k — average 
inter-dust grain spacing/dusty plasma Debye radius, indeed reveal dust solid face-centered cubic (fee), dust solid 
body-centered cubic (bec) and dust fluid phases for a set of T s and k values, as given in Ivlev et al. |16j where an 
empirical scaling for dust crystal melting is also given. 

The formation of dust Coulomb crystals and ordered dust structures have been observed in the sheath region of 
many laboratory experiments [17-20], where charged dust grains are kept together due to confining electrostatic 



2 



potentials in a plasma sheath. However, robust ordered dust structures may also be formed due to attractive forces 
[2] between negative dust grains associated with ion focusing and ion wakefields [21] [22] in a dusty plasma sheath 
with streaming ions, shadowing forces due to collisions with ions [221 HI], as well as due to overlapping Debye spheres 
[2"5] and dipole-dipole interactions [2~M2~B] . The alignment of charged dust grains in an assembly due to the attractive 
force associated with ion focusing and ion wakefield effects has been experimentally observed |29j . Furthermore, 
the collective behavior of dusty plasmas involving an ensembles of charged dust grains was recognized through the 
prediction of the dust acoustic wave (DAW) by Shukla [5U] at the First Capri Workshop on Dusty plasmas in May 
of 1989, where he suggested the existence of the nonlinear DAW in the presence of Boltzmann distributed electrons 
and ions, and massive, charged dust particles. This idea was then worked out in the first paper [5T] on the DAW. It 
must be stressed that there does not exist a counterpart of the DAW in an electron-ion plasma without charged dust 
grains, since the DAW is supported by the dust particle inertia, and the restoring force comes from the pressures of the 
inertialess hot electron and ions. Thus, similar to the Alfven wave in a magnetized plasma, the DAW is of fundamental 
importance in laboratory and space plasmas physics. The DAW is usually excited by an ion streaming instability, 
and has a frequency much smaller than the dusty plasma frequency, extending into the infra-sonic frequency range. 
Low- frequency (of the order of 10 Hz) DA fluctuations were first observed in the experiment of Chu et al. |17j . and 
have since been observed in many laboratory experiments world-wide [2] [5] [T7J 15214341 FI5] . and also in the Earth's 
ionosphere [35] . 

Ichimaru et al. |36j further extended the theory of strong coupling and viscosity coefficients for a high-density one 
component electron plasma. Berkovsky [37] developed a generalized hydrodynamic model for plasmas with strongly 
coupled ions and degenerate electrons, and used it to investigate the linear properties of modified ion-acoustic waves. 
A similar theory was developed for strongly correlated dust grains in dusty plasmas by Kaw and Sen 38J, who 
presented a generalized viscoelastic hydrodynamic model for strongly correlated dust grains and investigated the 
linear properties of dust acoustic waves, especially the low-frequency longitudinal and transverse modes in a strongly 
coupled dusty plasma. The latter model has also been extended to the weakly nonlinear regime |39j to study the 
propagation of small amplitude nonlinear dust acoustic waves in a strongly coupled dusty plasma. 

However, recently a number of laboratory experiments [4DH44] have reported observations of nonlinear DAWs in 
the form of extremely large amplitude DA shocks fiOTH2"] and DA solitary pulses [4"5tfi5] at kinetic levels. Physically, 
the large amplitude DA shocks are formed when nonlinearities in plasmas balance the DAW dissipation caused by the 
dust fluid viscosity coming from dust grain correlations in strongly coupled dusty plasmas, while DA solitary pulses 
arise in the collisionless regime due to the balance between the harmonic generation nonlinearities and the DAW 
dispersion. To the best of our knowledge, there are no theories for arbitrary large amplitude nonlinear, nonstationary 
DA shocks and DA solitary pulses in dusty plasmas with dust correlations. It should be stressed that small amplitude 
theories for DA shocks and DA solitary pulses based on the Burgers [46 , Korteweg-de Vries (KdV) , and KdV-Burgers 
equations 39^ are not suitable for explaining observations [551 I4TH44] that report anomalously high (up to 40% and 
beyond) dust density compressions. A large amplitude theory of Eliasson and Shukla [27] for a collisionless dusty 
plasma explains well the DAW steepening and nonlinear wave speed [41] [42] , but is unable to predict the shock width 
observed in the experiments. 

In this paper, we present a fully nonlinear, non-stationary unified theory for arbitrary large amplitude DA shocks 
and DA solitary pulses in a dusty plasma, taking into account the effects of strong coupling between charged dust 
grains, the nonlinear polarization force acting on charged dust grains due to thermal ions that shield negative dust 
grains, collisions between charged dust grains and neutrals, dust correlations decay rate, the dust fluid shear and bulk 
viscosities, etc. This gives a more complete picture of various non-ideal effects in dusty plasmas, and we are thus 
able to provide a comparison between our new non-stationary and fully nonlinear theory with the recent laboratory 
observations of DA shocks and DA solitary pulses [41TI44"] . Neglected are effects due to attractive forces (ion focusing, 
wake fields, etc.) between dust grains, which may affect the equation of state and transport coefficients of the system. 
These effects, however, are either small or depend on the moment transfer of streaming ions, which we do not consider 
here. 



II. MATHEMATICAL MODEL 

We consider a dusty plasma composed of inertialess electrons and ions, as well as strongly correlated negatively 
charged micron-sized dust particles of uniform sizes. In the presence of large amplitude ultra-low frequency DA waves, 
with uj <C v en , Vi n <C k 2 V£ e Til 10 ' where w is the wave frequency, v en (i^ n ) the electron (ion)-neutral collision frequency, 
k the wave number, and Vr e (Vri) the electron (ion) thermal speed. Both electrons and ions follow the Boltzmann 
law, since they can be considered inertialess on the timescale of the DAW period, and henceforth rapidly thermalize 
under the action of collisions. Thus, the electron and ion number densities are, respectively, n e = n e Q exp(e<f>/kBT e ), 
and rii = riio exp(— e^/fesTi), where n e o and are the unperturbed electron and ion number densities, respectively, 
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e the magnitude of the electron charge, <fi the electrostatic potential, ks the Boltzmann constant, and T e (Tj) the 
electron (ion) temperature. At equilibrium, we have the quasi-neutrality condition n^o = n e o + Zdtido, where Zd is 
the average number of electrons residing on a dust grain, and n d o the unperturbed dust number density. 

The dust particle dynamics associated with fully nonlinear, non-stationary DAWs in a strongly coupled dusty 
plasma is governed by the generalized hydrodynamic equations composed of the dust continuity equation (drid/dt) + 
V • {udYd) = 0, and the generalized dust momentum equation 
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taking into account finite amplitude convective and pressure nonlinearities |47j . nonlinear ion polarization force, strong 
dust coupling effects [36-381 4S], and dust neutral collisions [15] . Here d/dt = (dt/dt)+Vd-V is the total time derivative, 
n d and are the dust number density and dust fluid velocity, respectively, nid the dust mass, p d = nd,md the dust mass 
density, R = Zde 2 /AksTiXm is a parameter determining the effect of the polarization force [5D], which reduces the 
phase speed of the DAW, arising from interactions between thermal ions and negative dust grains, HdndksTd = P d the 
effective dust thermal pressure for a one component plasma (OCP) j3S], where pd = l + (l/3)u(T) + (T/9)du(r)/dT the 
compressibility, T = Z 2 e 2 / dksTd the ratio between the dust Coulomb and dust thermal energies, d = (S/inndo) 1 ^ 3 the 
Wigner-Seitz dust grain separation distance, and u(T) is a measure of the excess internal energy of the system, which 
reads [UGI2] «( r ) - -(V3/2)r 3 / 2 for T < 1 (viz. a liquid-like state), and u(T) = -O.SOr+O.QSr^+O.igr-^-O.Sl 
in a range 1 < T < 200. Furthermore, the effective dusty plasma Debye radius Ap = \De^Di/(^ e + ^Di) 1 ^ 2 , where 
Aue = (fcsT' e /47rn e oe 2 ) 1 / 2 and Xoi — (fcs7i/47rn i0 e 2 ) 1 / 2 are the ion and electron Debye radii, respectively. The 
dust-neutral collision frequency is given by the Epstein formula v dn = (8/3)^/2Trrn n n n r 2 VTn/''nd, where m n is 
the neutral mass, n n the neutral number density, the dust grain radius, Vr n — (kB^/uin) 1 / 2 the neutral thermal 
speed, and T n the neutral gas temperature. The visco-elastic properties of the dust fluids are characterized by the 
dust correlation relaxation time (351 [37] r r = [(£ + 4?7/3)/n ( 2o?d]/ [1 — Pd + 4u(r)/15], involving the shear and bulk 
viscosities rj and £, respectively. There are various approaches for calculating rj and £, which are widely discussed in 
the literature [52]. The DA wave potential <fi is obtained from Poisson's equation V 2 = 47re(n e — Uj + Z d rid), taking 
into account the dispersive effect due to the departure from quasi-neutrality. The ion drag force [5"3"H55] acting on 
a dust grain has been neglected in Eq. (1), which is justified since the ions are assumed to follow the Boltzmann 
distribution with no ion momentum flow. However, in a dusty plasmas with large dust particles and equilibrium ion 
flows, there can be an instability with a growth rate much smaller than the DAW frequency [2]. 



III. ONE-DIMENSIONAL QUASI-STATIONARY SHOCKS AND SOLITARY WAVES 

Let us now consider the simplest problem of one-dimensional nonlinear DAWs propagating along the x-axis in a 
Cartesian coordinate system. We define the dimensionless variables N — rid/rido, U = x • v^/Cd, and $ = e^/fc^Tj, 
where Cd = ^pd^D is the dust acoustic speed, uj p d = (AimdoZ^e 2 /rrido) 1 ^ 2 the dust plasma frequency, and x the unit 
vector along the a;— axis. We then have the dust continuity equation 
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where a — to pd T ri v — Vdn/^pd, D/DT = d/dT + Ud/dX, T = io p dt, X — x/\d, A = 
(£ + 4?7 /3) /m d ndoOJ P dd 2 (typical values [36 of (3 are roughly 1.04, 0.08, and 0.3 for T = 1, 10 and 160, respec- 
tively), Tq = [idTdl I Z d TiP , 7=1+ t(1 — P), P = Zdrido/nio, and t = Ti/T e . We are assuming here that the 
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constant parameter P is given for a set of experiments; however, it has been experimentally shown [5j that Zd is 
typically reduced for closely packed [d < Xd) dust grains. This effect, which can be important at high dust number 
densities, will be neglected here for simplicity. Furthermore, the dust charge fluctuation effect has been neglected, 
since the dust charging time-period (v~ ) is usually much shorter than the time period for the formation of nonlinear 
DAWs we are concerned with [5B], and the fugacity parameter F — 4irridoX Di rdV2 / + n e oTi/riioT e ) is smaller than 
1, where the expressions for v\ and are given in Refs. [2"1 [56], 

In a stationary frame such that all physical variables depend only on ( = X — MT with M = U/Cd, where U is 
the constant speed of the nonlinear DA waves, we have U — M(N — 1)/A, so that the dust momentum equation (3) 
reads 



! _ aM d 

+ A 9C 2 \NJ ' 
which couples with Poisson's equation 
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(5) 



a 2 $ 

7-^-2 = (1 - P) cxp(r$) - exp(-$) + PA. (6) 



Quasistationary DA shock waves exist only for v = 0, when the dust-neutral collisions can be neglected. Further- 
more, it is possible to derive a simple condition for the DA shock wave amplitudes depending on other parameters when 
the relaxation time for dust grain correlations is much smaller than the dust plasma period. Hence, for a = v = 0, 
Eq. (5) can be integrated once to obtain 

where we have used the boundary conditions N = 1 , $ = and d/d^ = at £ = +oo. The DA shock amplitude at 
C = — oo, where d/dt] = 0, N = N s i loc k > 1 and $ = & s hock < is now obtained from Eq. (7) as 

M 2 ( 1 \ iy 2'yR 

72 1 ) - ^ shock + -5-[l - exp(-$ shocfc /2)] + T lnA sflocfc = 0, (8) 
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while Eq. (6) yields 



exp(-$ shock ) - (1 - P) cxp(r$ s ^ oefc ) 
P 



(9) 



Using Eq. (9) we can eliminate N s f- l0 ck from Eq. (8) to obtain M as a function of the shock wave potential $> s hock 
for the parameters R, To, P, and r. The term proportional to /3/A in Eq. (7) works to smoothen the shock front, 
but does not influence the shock amplitude. The DA shocks are associated with a positive jump of the dust number 
density, N shock > 1, and a decrease of the potential, & s hock < 0, for M > C a , where C a = (f — R + To) 1 / 2 is the 
linear DAW speed in the long- wave limit d/d( = 0. Hence, the DA shocks are propagating with super- DA speeds in 
comparison with the upstream plasma. 



IV. COMPARISON WITH EXPERIMENTS 



In Figs. [T] and [2j we have used the plasma parameters of Refs. [4TJ [42] to study the nonlinear dynamics and the 
formation of shocks involving large amplitude DA pulses. The parameters of the experiment |4Tj are rij = 2 x I0 14 m~ 3 , 
Ti = 0.03eV, T e = 2.5eV, Z d = 2 x 10 3 , n d = 3 x 10 10 m- 3 , m d = ICT 15 kg, r d = 0.5 ptm, giving u> pd = 590s- 1 , 
Xd ~ 85/im, = 50mm/s, and d»2x 10~ 4 m. The used gas (argon, m n = 3.6 x I0 -29 kg) at the pressure 13 
Pa and temperature T n — 0.03 eV gives a neutral number density n n = 3 x 10 21 m~ 3 , and a dust-neutral collision 
frequency Vd n ~ 1 sT 1 ■ It was observed in the experiment [41] that a large amplitude dust density pulse self-steepened 
and formed a shock-like structure, which propagated with a mean speed of about 75mm/s, somewhat higher than the 
estimated dust acoustic speed. For the given parameters, we have A = 0.18, R = 0.28, P = 0.3, and r = 0.012. The 
normalized dust-neutral collision frequency v m 3 x 10 -3 is quite small, while the dust fluid viscosity due to strong 
dust coupling effects is more prominent. We choose /3 = 0.15, which is compatible with the experimental T > 1. 
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FIG. 1: The DA shock potential and associated dust number density as a function of M for P = 0.3, r = 0.012, R = 0.28, and 
To = 0.01. The DA shock potential is negative for increasing dust number density. The amplitudes increase with the increase 
of M. 




FIG. 2: The time and space evolution of (a) the dust number density and (b) the DA wave potential for a = 0.01, ft = 0.15, 
A = 0.18, v = 0.002, P = 0.3, R = 0.28, T = 0.01, and r = 0.012, corresponding to the plasma parameters of Refs. [4T1 142] . 



In addition, we choose a = T = 0.01. Figure [T] displays M as a function of dust number density and associated 
potential, obtained from Eqs. (8) and (9). In the small amplitude limit, viz. N s h oc k ~> 1 and $ shock — > 0, we have 
M — > C a ss 0.85. The DA shock speed M increases with increasing DA shock wave amplitudes, with an increase of the 
dust density and an associated negative potential. Figure [2] shows a simulation of the time-dependent system of Eqs. 
(5)-(8). As initial conditions, we used N = 1 + exp[-(A - 20) 2 /100] and U = 0.7exp[-(A - 20) 2 /100]. The profiles 



of the dust number density and DAW potential in Fig 
DA shock is formed, similar to the one in Fig. 5 of Ref. 



2^ show that the initial DA pulse steepens and a monotonic 
41] . The large amplitude (100%) dust density perturbations 
are associated with a negative potential <E ) « —0.25. The average speed of the DA density pulse is M w 1.4, in good 
agreement with Fig. 1 for N s h oc k ~ 2, M s» 1.3, and <fi s hock ~ —0.25. We found that monotonic (oscillatory) DA 
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shocks exist for f3 > A (f3 < A), and solitary waves in the limit j3 <C A. In dimensional units, the simulated nonlinear 
wave speed is about 70mm/s, which is close to the experimental value in Ref. |41j . 
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FIG. 3: a) The time and space evolution of the dust number density for a — f3 = v — R — To = 0, P = 0.51, and r = 0.025. b) 
The time variation of iV at X — 0. The driven DAW develops into spiky solitary DAW structures similar to those observed by 
Teng et al. [H]. 

We next turn to laboratory observations of large-amplitude localized DA solitary pulses in weakly collisional plasma 
discharges. Teng et al. jU] and Chang et al. [35] observed the formation of large amplitude, localized dust density 
structures, driven by a flow of ions towards the bottom of the plasma discharge. From the given parameters (44j [45] 
n e = 10 9 cm~ 3 , T e = 4eV, n t = 1.2 x 10 9 c m - 3 , T l = 0.05eV, Z d = 5000, n d w 3.7 x 10 4 cuT 3 (inter-dust distance 
about 0.3mm), and m,d = 6.9 x 10~ n g, we have uj p d = 200s -1 , Ad « 45 /jm, and Cd — 9mm/s. The observed 
nonlinear DA solitary pulses in Fig. 1(c) of Ref. jUj had a periodicity of about 2 mm, a mean speed of about 
45mm/s, and a crest width (measured at the height of N where N — 1) in the range 0.4 - 0.5, mm, with higher 
amplitude pulses having smaller widths. We believe that there are some uncertainties in the plasma parameters that 
could explain the relatively low value of Cd compared to the observed wave speed: Increasing the values of 2$ to 
0.10 eV and using the dust charging equation [e.g. Eq. (11) of Shukla and Eliasson [B]], we obtain Zd — 13800 for 
rii = 10 9 cm~ 3 and n e — 4.9 x 10 s cm -3 , giving uj p d = 540 s _1 , Ad = 76 /xm, and Cd — 41 mm/s, which is compatible 
with the experiment. Using these parameters in our model, we have P = 0.51 and r = 0.025, which we use in the 
simulation of the time-dependent system of equations (2)-(4). The results are displayed in Fig. [3] We drive the DAW 
resonantly by an external force of the form F = -0.01 sin[27r(A - T)/L] - 0.001 sin [2tt(X - T)/5L], added to the 
terms in the square parentheses in Eq. (3), where L = 26.3 is the observed wave periodicity (2 mm) normalized by 
Ad- The result in Fig. [3] shows almost periodic wave-trains that develop into narrow peaks, very similar to the ones 
observed by Teng et al. [44] . with density maxima about twice the ambient density and a typical width of about 4-5 
Debye radii corresponding to about 0.3-0.4 mm. These spikes may be interpreted as driven large amplitude solitary 
DAW structures due to a balance between the harmonic generation nonlinearities of the system and the dispersion 
provided by the departure from the quasi-neutrality condition. 

Bandyopadhyay et al. |43j studied how the speeds of DA solitary pulses depend on their amplitudes. The experimen- 
tal plasma parameters were n,- = 7 x 10 13 m~ 3 , T. t = 0.3 eV, T e = 8eV, n d = 10 10 m~ 3 , Z d = 3 X 10 3 , m d = 10~ 13 kg, 
giving uj p d = 51s -1 , A^> = 490/im and Cd — 25 mm/s, which corresponds to P = 0.43 and r = 0.038 in our model. 
The DA solitary pulses propagated with super-dust acoustic speeds, increasing with increasing amplitudes. In Fig. 
1(b) of Ref. [23] a pulse of 100% density amplitude propagates about 8 x 10~ 3 m in 0.24 s, giving a mean speed 
of Vd ~ 0.033 m/s, which corresponds to M = Vd/Cd — 1.33. Figure El shows a simulation result, where the initial 
condition consists of a wide pulse of the form N = l + 0.5exp[— (A — 20) /100], U = 0.5exp[-(A-20) 2 /100]. The DA 
pulse breaks up into three DA solitary wave structures propagating with the super-dust acoustic speed M > C a = 1. 
Small but finite amplitude DA solitary pulses have the density profile N = 1 + No sech 2 (C^ 2 ^/2), and the associated 
DAW potential $ = -(Af 2 PA / 7 )sech 2 (C 1/2 C/2}, where N = 3C ~//2BA1 2 P is the amplitude, C = l- 1/M 2 , and 
B = (I/27) [(1 — P)t 2 — 1 + 37 2 /M 4 P]. Figure Hb) exhibits that the numerically obtained amplitudes of the three 
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FIG. 4: a) The time and space evolution of the dust number density N for a — fi — u = R — To — 0, P = 0.43, and r = 0.038. 
The initial broad pulse breaks up into three separate DA solitary pulses propagating with the super-acoustic speed, similar to 
those observed by Bandyopadhyay et al. [43]. b) A comparison between the soliton amplitude obtained numerically (circles) 
with the theoretical amplitude No (solid line). 

DA solitary pulses compare favorably well with the theoretical amplitude AT . A pulse with a density amplitude of 
A = 2 would have a speed of M ~ 1.3, which is also in good agreement with the experiment of Ref. [35] . 

V. SUMMARY AND CONCLUSIONS 

In summary, we have presented a fully nonlinear, non-stationary unified theory for arbitrary large amplitudes DA 
shocks and DA solitary pulses in a strongly coupled dusty plasma. Our nonlinear theory is based on the Boltzmann 
distributed inertialess warm electrons and ions, Poisson's equation, the dust continuity equation, and the generalized 
viscoelastic dust momentum equation for strongly correlated charged dust grains. The governing nonlinear equations 
have been numerically solved to obtain the profiles of nonlinear DA waves, including the development of the DA shocks 
and DA solitary pulses. A comparison between our simulation results and recent experimental observations |411 142j 
of the DA shocks in laboratory dusty plasma discharges reveals a very good agreement with respect to the nonlinear 
DA wave speeds and DA shock wave smoothing due to strong coupling effects between charged dust particles. From 
the width of the DA shocks, one may, as suggested by Heinrich et al. 0T], infer the dust fluid viscosity. Furthermore, 
our simulation results of large amplitude DA solitary pulses also compare favorably well with the observations of 
Bandyopadhyay et al. (43] and Teng et al. [44]. Future experiments of nonlinear DAWs with higher precision 
measurements of the plasma parameters would be very valuable to benchmark the theoretical model. In closing, we 
stress that our fully nonlinear unified theory for DA shocks and DA solitary pulses remain valid for a dusty plasma 
with a weak magnetic field (of the order of 100 Gauss), since the latter is unable to magnetize micron-sized charged 
dust particles and would not affect the trajectories of electrons and ions that follow the Boltzmann law on the spatio- 
temporal scales of our interest. A weak magnetic field just provides confinement for the electrons, which are coupled 
with ions and negative dust grains through the space charge electric field of the DAW. 
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